Dynamical projection of atoms to Feshbach molecules at strong coupling 
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The dynamical atom/molecule projection, recently used to probe fermion pairing, is fast com- 
pared to collective fermion times, but slow on the Feshbach resonance width scale. Theory of 
detuning-induced dynamics of molecules coupled to resonantly associating atom pairs, employing 
a time-dependent many-body Green's function approach, is presented. An exact solution is found, 
predicting a 1/3 power law for molecule production efficiency at fast sweep. The results for s- and 
p-wave resonances are obtained and compared. The predicted production efficiency agrees with 
experimental observations for both condensed and incoherent molecules away from saturation. 
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Cold atomic Fermi gases, magnetically tuned to a Fes- 
hbach resonance region Q , host an intriguing strongly in- 
teracting many-body system Recently, pairing 
of fermions near the resonance was probed with the help 
of dynamical projection of atomic state on the molecular 
state 0,0 J achieved by a sweep through the resonance, 
followed by the detection of molecular Bose- Einstein con- 
densate. The sweep could be made very fast compared 
to typical fermion time scales, such as the collision fre- 
quency or inverse Fermi bandwidth and pairing energy 
gap, making the process a "snapshot probe" with regard 
to the collective fermion processes. 

On a single particle level, however, broad Feshbach 
resonances studied in Refs. 0, 0], exhibit strong atom- 
molecule coupling in a relatively wide detuning range. In 
this sense, the sweep speed 0, corresponds to essen- 
tially adiabatic atom/molecule conversion, slow on the 
scale of the resonance width. For example, in the JILA 
experiment , the Feshbach resonance width A_B ~ 10 G 
translates into Av ~ 180 MHz in detuning frequency, 
while the characteristic time of the magnetic field sweep, 
T, = {dt/dvY^^ ~ 1 /iS, is about 10^ times longer than 
Aiy^^. A similar estimate applies to the MIT fast pro- 
jection experiment Somewhat paradoxically, the 
fermions involved in this "slow" molecule formation are 
the same whose many-body state is being analyzed by the 
dynamical "snapshot" projection. Thus a correct physi- 
cal picture of the molecular state swept through the reso- 
nance must combine the adiabatic single-particle and the 
"snapshot" many-body aspects in a seamless way. 

Continuing efforts to use atom/molecule projection as 
investigative tool call for better understanding of the 
driven molecular state. The Landau-Zener model 
which fits the data well near saturation fl^, focuses on 
the adiabatic aspects, ignoring molecule dissociation into 
continuous spectrum of atom pairs. The dynamical mean 
field approach [HElElElIllEillSEil, which can be 
justified for bosons in the atomic BEC regime, lacks firm 
foundation in the fermion case. Recently, the many-body 
state overlap models 0, |^ were put forward. 
While providing some guidance, these approaches do not 
account for the experimentally relevant situation of broad 
resonance 0,0 when the "snapshot" many-body projec- 



tion is slow on the scale of individual molecule formation. 

Our objective is to describe molecules at a sweep 
fast compared to the elastic collisions, when only the 
quantum-mechanical processes are relevant. We develop 
a theoretical framework which accounts for resonance dis- 
sociation/association in the presence of time-dependent 
detuning as well as for fermion pairing correlations. We 
describe the molecules swept through the resonance us- 
ing a time-dependent Green's function which fully ac- 
counts free relative motion of the atoms associating to 
form molecules. While our method is quite general and 
applicable to Feshbach resonances with any angular mo- 
mentum, here we focus, for the sake of concreteness, on 
the s-wave case. We consider the evolution from equilib- 
rium at V = Uq, followed by an abrupt linear sweep: 
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with a the sweep rate. The g 


resonances 23, 24, 25| 
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;eneralization to the p-wave 
will be straightforward (see 
below) . 

Finding the time-dependent Green's function is a non- 
trivial mathematical problem, here solved exactly using 
an idea similar to that of the Wiener-Hopf method. The 
important time scale, characterizing the adiabaticity of 
the sweep (see Fig.Q] inset), is found to be 
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with g the atom- molecule coupling (see Eq.ljHl), and 
m the atom mass. The time scale tq can also be in- 
ferred, as noted by Altman and Vishwanath 22], from 



the adiabaticity condition lu 



< . 



for the time-dependent 



molecule energy huj. Different regimes arise depending 
on the relation between tq and vo/a, the time it takes 
the sweep to reach the resonance (Fig.^. The atom-to- 
molecule transformation takes place at times less than tq 
after crossing the resonance, where the evolution is nona- 
diabatic. At later times, the molecules, dressed by atom 
pairs, evolve adiabatically. For a fast sweep, aro 3> fo, 
the number of produced molecules scales with the sweep 
rate as a~^^^, while for slower sweep, aro <C I'o, the 
number of molecules scales as a~^. 
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Dyson equation '28'|, has the form 



Sweep rate, (^''/a)"^ 



1/3 



FIG. 1: Molecule number vs. the sweep rate (^X^/a 
at different initial detuning vq. The asymptotic regimes, 
Nm oc ct~^, a~^^^, correspond to slow and fast sweep. Inset: 
Molecular energy time dependence (0 with the time interval 
corresponding to nonadiabatic evolution marked. 



These results agree with the molecular number and 
condensate production efficiency reported by JILAgroup 
0. The sweep speeds \dt/dB\ « 10 - 80^s/G 's] cor- 
respond to vo/ce ~ 1 — 100 ^s with 1^0 — 0.1 — IG 
in the magnetic field units. The characteristic atom- 
molecule coupling « 1 GHz gives the adiabaticity time 
To ~ 10 — 20 /IS depending on the sweep speed. Thus 
with 0.2 < aTo/iyQ < 10 both the fast and the slow 
regimes are realized. Indeed, the molecule number ob- 
tained for different sweep speeds below saturation (see 
Fig. 5 in Ref. displaying the data for = 0.12 G) can 
be fitted quite accurately with the 1/3 power law depen- 
dence, Njn (x \dt/dB\^^^ , in agreement with our results. 
Also reasonable, by the order of magnitude, is the pre- 
dicted total number of produced molecules. Our conclu- 
sions regarding the incoherent molecule production chan- 
nel are consistent with the observed independence of the 
condensate fraction of the sweep speed. We obtain the 
same production efficiency for condensed and incoherent 
molecules fEa. (|22|l l. except near saturation. 

Let us recall the form of the two-channel Hamiltonian 
, describing pairs of fermions binding to form molecules 
at the resonance: 



H = Ha+ Ti-r, 



E 

p.p' 



(3) 



with Ha = J2pa 2^aito-ap^: = Ek(^ + ^)^k^k, and 
o-pa, the atom and molecule operators, a the spin {h = 
1). The detuning v is determined by the molecule and 
two-atom Zeeman energy difference, v = A/j, (B — Bq). 
The single molecule Green's function, obtained from 



G{LU,k) 



[b = uj- h iO , (4) 

4™ 



where = A(— w)^/^ is the self-energy describing 

molecule dissociation (s-wave), which arises after inte- 
grating over the 3d density of atom pair states A'^(e) oc 



1.' 



e^/^ along with a suitable ultraviolet regularization 

For time-independent v, the molecular state dressed 
by atom pairs, is described by the Green's function pole: 



Z{io) 



bj — u!{k) + iO 



(5) 



with Lo^k) given by (D — — v. Near the resonance, 

at <C AE^ — A^, neglecting lu compared to one 
obtains molecular energy quadratic in detuning: 



uj{k) = -{u/Xf + k^/Am. 



(6) 



At < 0, Eq.® gives the energy of molecules, while at 
> it describes a resonance in the two-fermion scat- 
tering mediated by virtual molecules l29l| . The residue Z 
defines the bare molecule weight in the physical molecule 
state, Z^^{uj) = dG^^/duj = l + |(-a))"i/2, which varies 
from zero to one across the resonance, at —AE^ ^ v < 0. 
At relatively small detuning, |i//Ai?»| 1, Z increases 
linearly: Z{uj) w 2\v\/X'^. 

To investigate molecule formation at the resonance, we 
consider the Green's function for the problem with time- 
dependent detuning i/{t). In this case, due to nonlocal 
character of S in the time domain, the molecule evolution 
is described by an integral-differential equation 17, 18] 

(idt - u{t) - |i) hk{t)-ji:k{t,t')hu{t')dt' = 77fe(t) (7) 

with7]fe(t) —g Je~*'^^VT (2;, a; the pairing field, 

and ^pa{x,t) — X]p '^p.o-e'^'^""''*. Here the self-energy is 

E.„>0^/..(.)e— '£ = f^, (8) 

and T,{t < t') vanishes due to the causality. 



The pairing field r/, which acts as a source in Eq.( |7| l. 
should be taken as a c-number for the condensed 
molecules (with fc = 0), and as an operator for the in- 
coherent molecules. Generally, its correlation function 
includes both the coherent and incoherent parts: 

{m,.m,<.) = [2TTf\m?5{u; - n)5{k) + k{lu, k), (9) 

where 770 describes a finite amplitude for two fermions 
to have opposite momenta in the paired state, with 
/i ^ 2Ep the chemical potential of a pair, and 
K{uj,k) — {{fjk,uj'>lk,u)) the dynamical pair correlator 
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which is nonzero even for ideal Fermi gas. We first con- 
sider the coherent molecule production, treating both 
V — f^oc"*^* and b{t) as c-numbers. The incoherent pair 
source K{uj, k) will be discussed subsequently below. 

The evolution problem |(7J) is non-elementary due to 
nonlocality of E(t,t'). Our approach employs an idea 
similar to that used in the Wiener-Hopf method. We first 
handle an auxiliary problem in which the linear sweep 
v{t) = vq — at extends from — oo to oo, and then modify 
the solution to describe the situation of interest 

The auxiliary problem in question is to find b(t), — oo < 
t < oo, which obeys a linear integral-differential equation 

{i2j-iyo~ T,{i2j) + at) b{t) = T]{t) , lo ^ idt (10) 

with a source term ri(t) of a generic form. It is conve- 
nient to go to Fourier representation, in which t = —id^j 
and the problem is reduced to an ordinary differential 
equation (w — S(w) — vq — iad^) b{uj) = rj(uj) for b(uj) = 
J e'^'^*b{t)dt. Solution of this equation, first order in d^, is 
found using the gauge transformation b{uj) e^'^'^'^^bluj) 
with the phase Lp satisfying 



— aLp' {ijj) = Dq{u}) — uj — S((jj) 
This problem is solved by the function 



Va- 



in) 



b{co) 



-ia-^e"^^'^^ 



e-"f^'^">T]{Lj')dLu'. (12) 



To verify l|12|l . one can compare its behavior to that 
of b{uj) = Ce*'^^'^-', the solution to the homogeneous 
problem H1Q(I . For uj large and positive, since {—lu — 
iO)^/^ — —iyJZj, we obtain the asymptotic behavior 
Do « if » -z(2A/3a)u;3/2, e±'^ oc e^""^"" 

(a = 2A/3a). Thus C = is required to eliminate expo- 
nential growth. Indeed, the asymptotic behavior of the 
integral in H12|l at large positive oj is non-exponential: 
/+°°e-*^('^')r;(w')rfw' « e-*v'('^)?7(w)/(/?'(w). (For any 
physical source, r]{LO oo) is algebraic.) At the same 
time, the behavior at large negative lo does not require 
special attention: </? is real at a; < 0, and so the expo- 
nentials e^"^ oscillate as e='=*'^(~'^)^^ without giving rise 
to "dangerous" asymptotic behavior. 

Now, having found the solution for the sweep span- 
ning the entire range — oo < i < -|-oo, let us consider 
the sweep trajectory (^. In this case, it is convenient to 
represent the function b{t) as a sum 6<(i) -I- b^(t), with 
6>_<(i) nonzero only at i > (i < 0), respectively, ob- 
tained by restricting b{t) on the half- line t > {t < 0). 
Then the evolution equation, in operator form written as 
(w — i/(t) — S(il')) b{t) = ri{t), can be represented as 



Dob< + {Do + at)by = ?/< -I- ??> 



(13) 



with 77>,< = 9(±t)ri[t) having the same meaning as 
6>,<(i), and Dq = Do{Ld) defined in Ea. ((TT|l . 



Let us project the terms on the left hand side on the re- 
gions t > 0, t < 0, taking into account the constraints due 
to causality. The integral operator E acts only forward, 
not backward in time, the property explicit in Eq.®. 
Thus the function (^Dq^uj) + at^ 6> is nonzero only at 
t > 0, while the function Do(t2;)5< has both the t > 
and the t < parts. This observation allows to write the 
problem as two separate problems for 6>,<(i) as follows: 



Dob 



< =V< 
< 



{Do + ai>; 



Dob 



oo< =V>, 
> 



(14) 



where [...]<,> denotes the part of the function at i > 
(t < 0), with zero value on the opposite half- line. Now, 
we can solve the first equation for 6< and substitute the 
result in the second equation, which (after some algebra) 
can be brought to the form 



(Do + at)by = Do 



Do'v 



(15) 



We note that 6> and the function on the right-hand side 
are both nonzero only at t > 0. This allows to treat this 
equation as Ea. l|lU|l . formally extending the linear time 
dependence at to negative t. Using the above result, we 
obtain the answer in Fourier representation of the form 
(|12|l with 77 replaced by 



fi{uj) = Do bo^il 



Do{io) 



D^^{uj')ri{uj')dw' 
6 — i{uj — uj') 27r 



Now, let us consider the source r]{t) = jyoe^*^*, describ- 
ing coherent fermion pairs with the chemical potential 
/i/2 per particle. In this case, 'q{uj) = 2'K7]o5{lo — fi) 
and fi{uj) = Dq^ {iJ,)rioDo{uj) / {S — i{uj — fi)). Inserting 
f] in Ea. (|12|) . and using the identity (6 — i(uj ^ ^))^^ = 
j°° gi(w-^)T^^^ gj^j^ ^ closed form representation 

bH ^ , "l"^" , + Ab(c.), 



S — i{uj — n) 



(16) 



-t-00 



with A = Dq^(ii). (To obtain l|16(l . we transformed 
the integral over uj' by writing Do{lj') — ad{Lj'T — 
ip{oj')) / duj' — ar and integrating by parts.) Since the 
first term of lfTB|) gives the would be b{u!) in the absence 
of the sweep, A6(w) describes the effect of the sweep. 

Now, let us analyze the asymptotic behavior of b{t) ~ 
J e~*"*&(a;)dw/27r at large positive t :§> tq, lyo/a. In this 
case, the integral over lo is controlled by large negative uj, 
which can be seen with the help of the stationary phase 
approximation. Indeed, the saddle point of ~ujt + 
(p{uj), obtained from ip' = t, ait ^ +00 implies uj — > —00. 
With that in mind, we obtain the asymptotic for b{t) by 
setting the lower integration limit in Ea. H16|) at uj — —00, 
leading to the central result of this work: 



Ab{t) 



TF{T)dT, (17) 
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+ CXD 



(18) 



The quahtative behavior of F(t) can be analyzed using 
the stationary phase approximation. We obtain F(t) = 
•\/— 27ri/(^"(a;*)e*"*'^~''^^'^*\ where the stationary phase 
equation for w*, given by Dq{uj^:) + at = 0, has a real so- 
lution = — a^(t — to)^/A^ only for t > to = I'o/a. Re- 
lating the curvature tp" and the Green's function residue, 
yields the asymptotic form 



-aip" = D'o = Z~ 



F{t > to) = (27rmZ(w,)) 



1/2, 



(19) 



with Z(cj,) = 2a(t-to)/A2. (The self-energy -dominated 
-Do(a;) — — z/Q — S(ijj), appropriate for broad s-wave Fesh- 
bach resonance, was used in the above estimates.) Thus 
F(t) grows as (t — to)^/^ and oscillates at (t — to)/To » 1, 
decreasing exponentially at t — to < 0. 

To apply Eq. H17|l to the experimental situation we take 
into account that /i ^ T^/tq, H/tQ. (Indeed, /i ^ 2Ep, 
with Ep = 0.35 = 50 KHz in Ref. 0.) We evaluate 
the integral in H17|) using the stationary phase form (|19|l : 

Afe(t) - (2aZ(t))i/2A77oe'^(*-*°) Vi^o + C2T0) , (20) 

Ci = (i/3)i/2r(l/2), C2 = (i/3)i/^r(5/6). The asymp- 
totic number of molecules at time t is evaluated as TVm ^ = 
|A&(t)p (see Fig0. The fast and slow sweep regimes can 
be identified: at to <C tq we obtain Nm oc a~^^^, while 
at to ^ To we have Nm cc a~^. 

The incoherent molecule production can be studied in 
a similar manner. Using the operator rjk as a source in 
© , and averaging over its dynamical correlations we 
obtain the molecule momentum distribution 

oc 

^™(fc)=^|^^^*(t) Je-'^^TF{T)dT\\{cu,k). (21) 



with A = Dq^{lj). As a function of frequency, K is 
nonzero at a; < 2Ep. At Ep ^ S/to, ?i/to, the expres- 
sion is w-independent, as above. Factoring it out, 
we conclude that the condensed and incoherent molecule 
production efficiencies are identical. The molecule con- 
densate fraction is then expressed through the fermion 
pair fraction: 



N, 



(0) 



(22) 



fia-ix) — •4)cr(x)'ij}a{x). We note that the incoherent con- 
tribution exists even in the absence of pairing. For ideal 
fermions at density n, we have K(uj, k) = ^g'^n{l — 
u)^{l + u/2)9{l — u),u = k/2pp, which corresponds to a 
broad molecule momentum distribution with k < 2pp. 

The approach presented above yields accurate results 
for the atom/molecule projection in a wide range of sweep 



rates, fast and slow, as long as the times To, to are short 
on the scale of Ep. The only limitation stems from the 
assumption of a steady source, which describes the situa- 
tion when the fraction of atoms converted into molecules 
is small. The depletion effects, which are different for the 
condensed and incoherent molecules, can be incorporated 
in the framework of a quantum kinetic equation. 

The above method is applicable to the p-wave Feshbach 
resonance case, with the essential modification in the self- 
energy form S](t^) oc (-w)^/^ j^^g^. 

This self-energy 

is an irrelevant perturbation near the resonance, so we 
have the atom/molecule conversion time To = a~^^^ and 
Z = 1, as for weak coupling. Thus one can use Do{uj) = 
Lu — fo in limi , yielding the result identical to H17(l with 



F(t) = 



.(t-to)W/2a^^ = (27rza)i/2e-»"(*-*o)V2_ 



The number of produced molecules, both condensed and 
incoherent, scales as inverse sweep rate a^^. The pro- 
duction is less efficient than in the s-channel case due to 
weaker coupling at resonance. 

In summary, molecule production at Feshbach reso- 
nance is considered as a many-body problem for which 
the exact Green's function is obtained using Wiener-Hopf 
method. The theory is applied to the s-wave and p-wave 
resonances. The slow and fast sweep regimes are identi- 
fied in the s-wave case, controled by the adiabaticity time 
scale (O. The predicted power law 1/3 for the molecule 
production, as well as the total molecule number, are 
found to be in agreement with observations away from 
saturation 6]. The independence of the produced con- 
densate fraction on the sweep rate observed at fast sweep 
is also explained by this theory. 

We are grateful to Dmitry Petrov for useful comments. 
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